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Abstract 
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it without solving any differential equations. Examples of ODEs not having point symmetries are 
shown to be solvable using this algorithm. The scheme was implemented in Maple, in the framework 
of the ODEtools package and its ODE-solver. A comparison between this implementation and other 
computer algebra ODE-solvers in tackling non-linear examples from Kamke's book is shown. 
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1 Introduction 



^From a practical point of view, when developing solving methods for ODEs, what one actually does is 
attempt to determine families of ODEs which can be transformed into algebraic problems or into simple 
ODEs such as0 y' = F(x) or y' — F(y) by changes of variables or equivalent processes. For high order 
ODEs, one hopes that such a simplification of the problem will be possible after successive reductions of 
order. Some more powerful schemes are also able to exploit other information, as for instance integrating 
factors or the ODE's symmetries, and so to try a multiple reduction of order at once (see for instance [0 
and |). 

In the specific case of integrating factors, although in principle it can always be determined whether 
a given ODE is exact (a total derivative), there is no known universal scheme for making ODEs exact. 
Actually, for reorder ODEs - as in the case of symmetries - integrating factors are determined as solutions 
of an n t,l order linear PDE in n+1 variables, and to solve this determining PDE is a major problem in 
itself. 

Bearing this in mind, this paper presents a method for obtaining integrating factors of the form 
/j,(x, y') and fx{y, y') for non-linear second order explicit ODEs[], using a different approach, based only 
on a computerized analysis of the pattern of the given ODE. That is, for a given ODE, if an integrating 
factor with such a functional dependence exists, the scheme returns the integrating factor itself without 
solving any differential equations. 

The exposition is organized as follows. In sec. 0, the use of integrating factors for solving ODEs is 
briefly reviewed. In sec. ^, the scheme for obtaining the aforementioned integrating factors fi(x,y') or 
y') is presented and some examples are given. In sec. ||, some aspects of the integrating factor and 
symmetry approaches are reviewed, and their complementariness is illustrated with two ODE families 
not having point symmetries. Sec. || contains some statistics concerning the new solving method and 
the second order non-linear ODEs found in Kamke's book, as well as a comparison of performances of 
computer algebra packages in solving a related subset of these ODEs. In sec. || the computer algebra 
implementation of the scheme in the framework of the ODEtools package j| is outlined, and a description 
of the package's new command, redode, is presented. Finally, the conclusions contain some general 
remarks about this work and its possible extensions. 

Aside from this, in the Appendix, a table containing extra information concerning integrating factors 
for some of Kamke's ODEs is presented. 

2 Integrating factors and reductions of order 
2.1 First order ODEs 

The idea of looking for an integrating factor (/x) is usually presented in the framework of solving a given 
first order ODE, say, 

y' = $(x,y) (I) 
If by multiplying Eq.([l]) by a factor [i(x,y), the ODE becomes a total derivative^, 

»(x,y)(y'-*(x,y)) = -^R(x,y) (2) 
for some function R, then one can look for /i as a solution to the first order PDE: 

£ + = M-0 (3) 



dx dy 

which arises as the exactness condition for the problem (see Eq. ([?])). Once [i has been obtained, R(x, y) 
- an implicit form solution - can be calculated as a line integral. 



3 Throughout this article, we use the notation y = y(x),y' = ^Jri?/™' = ^ppr- 

4 We say that a second order ODE is in explicit form when it appears as y" — &(x, y, y') = 0. 

5 In this paper we use the term "integrating factor" in connection with the explicit form of the ODE, i.e., the ODE, 
turned exact by taking the product fiODE, is assumed to be of the form y" = &(x,y,y') or y" — &(x,y,y') = 0. 
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Although to solve Eq.([|) for \i is as difficult as the original problem, it turns out that for a given 

N(x,y), when a solution of the form /j,(x,y) = A* (q) N(x,y) exists - q is either x or y only - /z can 

be determined by solving an auxiliary linear first order ODE. For example, introducing fj,(x,y) = A* 
(x) N(x, y) and M(x, y) = N(x, y) &(x, y), one obtains: 

1 fdM dN\ 

H(x) = Cie V 7 (4) 

and a solution A 1 (x) exists only when the integrand in above does not depend on y. This gives both an 
existence condition and an explicit solution to the problem; however, the advantages of the scheme arc 
only apparent since there is no way to determine in advance what would be the appropriate N(x,y). 

2.2 High order ODEs 

Integrating factors for high order ODEs are defined as in the first order case. Here, we consider 
li(x, y, y' , j/" -1 **) to be an integrating factor for an n t/l order ODE, say 

yW = ^(x,y,y',...,y^) (5) 
if after multiplying the explicit ODE by // we obtain a total derivative: 

^ (n) -*)=S (6) 

for some function R(x, y, y' , j/™ -1 -*). To determine fi, one can try to solve for it in the exactness 
condition, obtained applying Euler's operator to the total derivative H = /i (j/ n ' — $) : 

+ (-1)"— ( —) = (7) 
Now, it can be shown by induction that Eq. (Q) is always of the form 

A(x, y, y', ...,y (2n ^) + y (2n ~ 2) B(x, y, y>, y^) = (8) 

where A is of degree n-1 in y( n > and linear in z/ fe ) for n < k < (2n — 3), so that Eq.(@) can be split into 
an overdetermined system of PDEs for fi. For example, for second order ODEs Eq.(|7]) is of the form 

A(x,y,y')+y" B(x,y,y') =0 (9) 
Hence, by taking A(x, y, y') = and B(x, y, y') = we have a system of two PDEs for /x: 

A(x,y,y')= (fO) 
d 2 u ( d 2 n \ , <9m\ x ( <9 2 $ <9$ / <9 2 $ \ A 



dH 




(-) 


d 2 


(—) 


dy 


dx ' 


[dy'J 


dx 2 


\dy") 



dy'dx \dy'dy J dy J \dy'dx dy \dy'dy 
dy J \\oy J dy \dy J dy' \dxdy J J 



d^\ 






<9$ d 2 [i 


dy 1 ) 




\dxj 


1 dy' dx 2 



= 



B(x,y,y')= (11) 
fd 2 fi \ ^ n fd(i\ d<S> 9 2 $ d 2 [i ( d 2 [i\ , „ 



dy \dy' 2 J \dy'J dy' r dy' 2 dy'dx \dy'dy 
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Nonetheless, there are no general rules which might help in solving these PDEs[]. 

Alternatively, a possible strategy for directly obtaining R instead of looking for /j, can be formulated 
as follows. Consider the first order linear operator associated to Eq.(||) 

where x, y and its derivatives are all treated as independent variables on the same footing. Now 

A(R) = (13) 

and there are n functionally independent solutions (first integrals) to the problem. In some cases, a first 
integral R such that dR/dy^ 1 ^ ^ can be obtained as the solution to a subset of the characteristic 
strip of A(R) = 0, or by other means. 

3 Integrating Factors and ODE patterns 

Since the classical way for determining integrating factors leads to problems similar in difficulty to solving 
the ODEs themselves, we consider here a different approach, based on a careful matching of an ODE 
pattern. 

The starting point is the observation that it is trivial to solve the inverse problem; i.e., to find the 
most general ODE having a given /i. In fact, from Eq.(^), we have 

»(x,y,y>,...,y^) = -^ (14) 
and hence the reduced ODE R is of the form 

R = G(x, y, y(»- 2 )) + f /x dy^ (15) 



for some function G. Inserting Eq.(|l5|) into Eq.(|l3|) and solving for leads to the general form of an 
ODE having /i as integrating factor: 



y (n) = _± 

A 6 



Yx (/ ^ (W_1) + G ) + - + (/ ^ (W_1) + G ) 



(16) 



The expression above then becomes an ODE pattern which one can try to match against an input ODE. 
The generation of such pattern matching routines is difficult, even for restricted subfamilies of integrating 
factor, but once built, they are a powerful and computationally efficient way to reduce the order of the 
corresponding ODEs (see sec. ||). 

3.1 Second order ODEs and the integrating factor family fi(x,y') 

In the case of second order ODEs, if instead of considering the general case fx(x, y, y') we restrict the 
family of integrating factors under consideration to fi(x,y'), Eq.([l5|) - the reduced ODE - becomes 

R(x,y,y')=F(x,y') + G(x,y) (17) 

for some functions G and F, where 

l i{x,y') = F yl {x,y') (18) 
(we denote Fyi — 7p7-)- Eq. (flit) can then be written in terms of F and G as 

6 In a recent work by H (1997), the authors arrive at Eq.(0) and Eq.(|llj) - numbered there as 13-5) and (3.8) - departing 
from the adjoint linearized system corresponding to a given ODE; the possible splitting of Eq. (pi) into an overdetermined 
system for /i is also mentioned. However, in formula (3.5) °^jjj' la * W< jj^j' °^^ c ''^ ^h nv^ annears^ replaced by "3?(a;, y, y'), 



and the authors discuss possible alternatives to tackle Eqs.(^) and ( |ll| ) instead of Eqs.(|10|) and (p.l|). 
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y" = *(x,y,y>) ee _F^v') + G x {xy) + G y (x,y)y> (lg) 

Fyi \ x ,y ) 

The idea is now to build a routine to determine if a given ODE can be written in the form Eq.(p^), in 
which case it will have an integrating factor of the form fj,(x,y'), and if so, determine it, leading to the 
reduced ODE Eq.jl^) by means of standard methods (see for instance || p. 221). The feasibility of such 
a computational routine is based on the following theorem. 

Theorem 1 Given a nonlinear second order ODE 

y" = $(x,y,y') (20) 

7^ 0) Q f or which an integrating factor of the form fi(x, y') exists, such an integrating factor can be 
systematically determined without solving any differential equations. 



PROOF. We divide the proof in two steps. In the first step we assume that, given Eq.(20), it is always 



possible to determine fj,(x 7 y') up to a factor depending on x; that is, to find some T(x,y') satisfying 

T M ) = ^ (2D 
fi(x) 

for some unknown function l^(x). We then prove that the knowledge of T(x,y') is enough to determine 

by means of a simple integral, hence leading to the desired fj>(x, y'). 
In a second step, we prove our assumption, that is, we show how to find T{x,y') satisfying Eq. (pi)) , 
concluding the proof of the theorem. 

3.1.1 Determination of fi>(x) when !F(x,y') is known 

Starting with the first aforementioned step, we assume that we can determine Tix^y 1 ). It follows from 
Eqs.©, © and @ that 

± ( *(», y, y') ?{x, y')) = °^ v) ± G ^ x > v) V - (22) 
so that by taking coefficients of y' in we obtain 



<fH = $ y {x,y,y')F{x,y')-y' ($y(x,y,y') F(x,y')J 

® ( ^ f„ ». „./\ -TV™ „./\ 1 _ Gyy{x,y) 



G yx (x, y) 
Mas) 

r2 = -^ 7 [<t> y (x,y,y')?(x, y >)) = ~™^ 1 " (23) 

°y v / n(x) 



Similarly, we obtain 



d (,, , s „, „\ F y > x ( x ,y') + G y( x >y) 



<P3 = -jr-i ($(x, y ,y') F(x,y')) = 

Q y \ ) ii{ x ) 

- &M)= (24) 
°V At (a:) 



Now, since Eq.(|20|) is nonlinear by hypothesis, either ip2 or (p4 is different from zero, so that at least 

one of the pairs of ratios {<pi, ^2} or {ip^, 1/74} can be used to determine A 1 (x) as the solution of an 
auxiliary first order linear ODE. For example, if tp 2 7^ 0, 

7 ODEs missing y may also have integrating factors of the form fj,(x,y'), which cannot be determined using the scheme 
here presented. However, such integrating factors are not really relevant since these ODEs can always be reduced to first 
order by a simple change of variables. 
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(tpi(x,y) = ^- (ifrfay) (25) 



and we obtain 



/ 



1 / <9</?i dtp 2 



\ ■ , i / rfa; 

, LDo \ OV OX I 

H{x)=e~ ^ 2 V ^ y (26) 



If y>2 = then (/34 7^ and we obtain 



1 / dips dip. 



, 094 V 9y' dx I 
=e ^ 4 V y 7 (27) 

Eqs.(p6|) and ( p7j ) alternatively give both an explicit solution to the problem and an existence condition, 

since a solution M (x) - and hence an integrating factor of the form n{x,y') - exists if the integrand in 
Eq.(p6|) or Eq.(p7j) only depends on x. A 

Example: Kamke's ODE 37. 

y" = -2yy' - f(x)(y' + y 2 )+g(x) (28) 
This example is interesting^] because it has no po int sy mmetries for arbitrary f(x) and g(x) (see sec. [|). 



For this ODE, !F(x,y') was determined (see sec. 3.1.2 ) as: 



F(x,y') = l (29) 

from which (Eq.©) 

G yx (x,y) + G yy (x,y)y' 



M(z) 



-2yf(x)-2y' (30) 



and then as in Eq.(p3[) we obtain 



<Pi = -2yf(x) 

ip 2 = -2 (31) 



Using this in Eq. fl26|) , we get 

/ }{x) dx 

V(x)=e J (32) 
and so, from Eq.(pl|), since !F(x,y') = 1, fj,(x,y') =fJ-(x). 

3.1.2 Determination of J 7 (x, y') 

We now prove our assumption, that is, we show how to obtain a function !F(x, y') satisfying Eq.(pl|) from 
the knowledge of <i>(x,y,y'), and without solving any differential equations. 

Since we have already assumed that the given ODE has an integrating factor of the form /i(x,y'), 
then there exist some functions F(x,y') and G(x,y) such that it is possible to rewrite $(x, y, y') - the 
right-hand-side (RHS) of the given ODE - as in Eq.([l9|). We then start by considering the expression 

^_d$_ G xy (x,y) + Gyy(x,y) y' 

~ dy ~ F yl (x,y') [d6) 



8 For ODE 6.37, Kamke shows a reduction of order_to a general Riccati ODE, based on the theory for ODEs having 
solutions with no movable critical points - sec W, and fel p 331. 



G 



and the possible cases. 
Case A 

The first case happens when the ratio G xv (x,y)/G yy (x,y) depends on y; i.e., G xy {x,y) and G yy (x,y) 
are linearly independent w.r.t y. To determine whether this is the case, note that we cannot just analyze 
the mentioned ratio itself since it is unknown. However, we can always select the factors of T containing 
y, and check if this expression also contains y' . If so, we just determine F i {x, y') up to a factor depending 

on x, that is, the required J-(x, y ), as the reciprocal of the factors of T which depend on y' but not y. A 
Example: Kamke's ODE 226 

This ODE is presented in Kamke's book already in exact form, so we start by rewriting it in explicit 
form as 

,, x 2 yy' + xy 2 .„,. 

y" = yy , y (34) 



y 



We determine T (Eq.p^)) as 



x(xy' + 2y) . „. 

T = v y - — ^ (35) 



The only factor of T containing y is: 



xy' + 2y (36) 
and since this also depends on y', T{x,y') is immediately given by 

T(x,y') = y> (37) 

Case B 

When the expression formed by all the factors of T containing y does not contain y' , it is impossible 
to determine a priori whether one of the functions {G xy (x, y), G yy (x,y)} is zero, or alternatively their 
ratio does not depend on y. We then proceed by assuming the former, build an expression for !F(x, y') 

as in Case A, and determine A* (x) as explained in the previous subsection. If this doesn't lead to the 
required integrating factor, we then proceed as follows. 

Case C 

In this case, we assume that neither G xy (x,y) nor G yy (x,y) are zero and their ratio is a function of 
just x, so that we have 



G xy (x,y) = v 1 (x) w(x,y) 

G yy (x,y) = v 2 (x) w(x, y) (38) 
for some unknown functions v\ (x) and v 2 (x) , such that Eq. ( |33| ) can be factored as 

r = M*,y) (vi{ i + ( V2 ^ y,) (39) 

for some function w(x,y) which can always be determined as the factors of T depending on y. To 
determine Fyi(x, y') up to a factor depending on x, we then need to determine the ratio vi (x)/«2 [%) ■ For 
this purpose, from Eq.(|38|) wc build an auxiliary PDE for G v (x,y), 

G xy (x,y) = ^G y , y (x,y) (40) 
The general solution of Eq.([40|) is given by 

G y (x,y)=G(y + p(x)) (41) 
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where Q is an arbitrary function of its argument and for convenience we introduced 

p'(x) = v 1 (x)/v 2 (x) (42) 

We can now determine p'(x), that is, the ratio vi/v 2 we were looking for, as follows. Taking into account 
Eq.(p8|), we arrive at 



v 2 (x) w(x,y) = Q'(y + p(x)) (43) 
By taking the ratio between this expression and its derivative w.r.t y we obtain 

-ui _l ( \\- dw I S"(V + P( X )) f a a\ 

n{y+p{ X )) = ^/ W = g>{y+p{x)) (44) 

that is, a function of y + p(x) only, which we can determine since we know w{x, y). If Ti.' ^ 0, we obtain 
p'{x) as 

/ d 2 w \ / dm \ dm 

_ mm _ W-Uw 
dx dy (0) »-(f) 

Once we determined p' (x), from Eq.(|39|) we determine J-(x, y') as 

H*,y') = ^f^ (46) 

where T, w(x,y) and p'{x) are now all known. A 
Example: Kamke's ODE 136. 

We begin by writing the ODE in explicit form as 

y" = (47) 

x — y 

This example is interesting since the standard search for point symmetries is frustrated from the very 
beginning: the determining PDE for the problem will not split due to the presence of an arbitrary function 
of y'. Here T (Eq.(|33"|)) is determined as 

r = - 7 ^ T2 (48) 

(x - yy 



and w(x, y) as 



Then H(y + p(x)) (Eq.@) becomes 



and hence, from Eq.fl45|), p'{x) is 



so from Eq.(46): 



Case D 



w{x,y) = 1 (49) 
(x - y) z 



H = — ^— (50) 
x-y 



p'{x) = -1 (51) 



^ = W) (52) 



We now discuss how to obtain p'(x) when 7t'(y + p(x)) = 0. We consider at first the case in which 
H = 0, hence Q" = 0, and so, recalling Eq.(41), we see that 



8 



G(x, y) =B 1 (y+ p(x)) 2 + B 2 (y+ p(x)) + g(x) (53) 
for some function g(x) and some constants B\, B 2 - Recalling Eq. (|l9|) , &(x,y,y') takes the form 

/\ _ v') + 9_M + ( 2i?1 ^ + + B 2 )(y' +P'(x)) 

v{x,y,y)- F^W) (54j 

We can now obtain explicit equations where the only unknown is p{x) as follows. First, from the knowledge 
of T and $ we build the two explicit expressions: 

1 F y' 

A ^T = ~2B 1 (yi+p>(x)) (55) 

and 

* s - y = J-y V> + p{x) + §■ (56) 

It is now clear from Eq.(^5|) and Eq.(p6|) that A and & are related by the following equation: 



d_ 

dx 



({y'+p'{x)) a) + A (V + + ^ (57) 

where the only unknowns are p{x), B\, and B 2 . By differentiating the equation above w.r.t y' and x we 
obtain two equations where the only unknown is p'(x): 



A y ,p"(x) + (A xy , + * y , y ,)(y'+p'(x)) + A X + 2«y = (58) 

A p"'(x) + {A xx + * y , x )(y' + p'(x)) + (A x + * y Ap"(x) + * x =p'{x) (59) 

As a shortcut, if (A / + ^ i A/A i depends on y', then we can build a linear algebraic equation for 
p'{x) by solving for p"(x) in Eq.(J5q) and differentiating w.r.t. y' . Otherwise, in general we obtain p'(x) 
by solving a linear algebraic equation built by eliminating p"{x) between Eq.([58|) and Eq. (|59])P|. 



If Eq.(58) depends neither on p'{x) nor on p"(x) this scheme will not succeed. However, it is possible 
to prove that in that case the original ODE is already linear, and easy to solve. To see this, we set to 
zero the coefficients of p'(x) and p"(x) in Eq.(|5^), obtaining: 

V = V + *2/V =A * + 2 *2/' = ° ( 6 °) 



from which A is a function of x only, and then 



*y>y> = (61) 



If we now rewrite F(x,y') as in 



F(x, y') = Z(x, y') - g{x) - A(y' + p'f B 1 (62) 

and introduce this expression in Eq.(|55|), we obtain Z i — 0; similarly, using this result, Eq.(|56]), Eq.(|6l|) 

and Eq.(|6^) we obtain Z x = 0. Hence, Z is a constant. Finally, taking into account that Z is constant, 
Eq.(|62|) and Eq.([5^), we see that the ODE Eq.(|2C|) which led us to this case is just a non-homogeneous 
linear ODE of the form 

(y + p)" + (A'(y + P )' - 2(y + p) - B 2 /B 1 ) /2A = (63) 
which homogeneous part does not depend on p(x): 

9 From Eq. (|)5|), A ^ 0, so that Eq. (|5S|) always depends on v"'(x). and solving Eq.(^fj) for p"(x) and substituting twice 
into Eq. will lead to the desired equation for p'(x). If Eq.(|58[) depends on p'(x) but not on p"(x), then Eq. (pi|) itself is 
already a linear algebraic equation for p'(x). 
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and which solution is in any case straightforward. A 
Example: Kamke's ODE 66. 
This ODE is given by 

y" = a{c + bx + y)(y' 2 + lf' 2 (65) 
Proceeding as in Case A, we determine T, w(x,y), and 7i{y + p(x)) as 

T = a (V 2 + 1) 3/2 ; w(x,y) = l; H = (66) 



From the last equation we realize that we are in Case D. We determine A and (Eqs. (pq), (p6[)) as: 

1 



A = 



, 3/2 

a 



We then build Eq.(p7f) for this ODE: 



* = c + bx (67) 



P ' l{ - X) +c + bx = p(x) + J^- (68) 



3 , - y [% 2 = (69) 



FM)= , ,3/2 ( 71 ) 



/ /2 i A 3 ^ 2 2S X 
(/ +1J a 

Differentiating w.r.t. y' leads to Eq.(|5g|): 

P"(*)y' 
(,' 2 + i) 5 " a 

from which it follows that p"(x) = 0. Using this in Eq.(p)9|) we obtain: 

p'(x) = b (70) 

after which Eq.([l6|) becomes 

v' + b 
a (y' 2 + l) 

Case E 

We now show how to obtain p'{x) when Ti'{y + p{x)) — and H = Q" /Q' is a constant, Bi, which is 
different from zero; so Q' is an exponential function of its argument (y + p(x)) and hence from Eq.(|4l|) 

G(x, y) = B 2 e^ +p(x))B ' + (y + p(x))B 3 + g(x) (72) 

for some constants B 2 , -B3 and some function g{x). In this case, it is always possible to arrive at an 
algebraic equation for p'(x), though the case entails some subtleties. First of all, $(x, y, y') will be of the 
form 

Af » _ F x (x,y')+g'(x) + (j? 2 j? ie fa+^»^ + B 3 ) (y'+ P '(x)) 

9(x,y, V) - F ^ {xyl) [16) 

Now, taking advantage of the fact that we explicitly know B\ , we build our first explicit expression by 
dividing B ie yBl by T: 
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A = 2- (74) 

B 2 ePW B i (y'+p'(x)) K ' 

We now multiply $ by A and subtract B\e BlV to obtain our second explicit expression: 



B 2 eP ( - x '> Bl V y' +p'{x) 
Now, as in Case D, A and *f> are related by 



d_ 

dx 



(y' + p'(x)) a) + (y> + p'(x)) p\x)AB 1 + A f(y' + p\x)) *j 



B ? eP( x ) B i 



(76) 



where the only unknowns are B 2 , B$ and p(x). We build a first equation for p'(x) by differentiating 
Eq.(f76|) with respect to y 1 



(p"(x) + P '(x) 2 B^j A y , +p'(x) (y'A^Bx + AB l + A xy , + ^ y , y , 



-2y + A x + y 'A xyl + y'y yly , = 



(77) 



The problem now is that, due to the exponential on the RHS of Eq.(|76|), differently from Case D, we are 
not able to obtain a second expression for p'(x) by differentiating w.r.t x. The alternative we have found 
to determine p'(x) can be summarized as follows. 

We first note that if Ay, = 0, Eq.(^) is already a linear algebraic equation for p'fj, so that we are 
only worried with the case Ay, ^ 0. With this in mind, we divide Eq.(^) by A y , and, if the resulting 
expression depends on y' ', we directly obtain a linear algebraic equation in p'{x) just differentiating w.r.t 
y'.A 

Example: 

y or + y — 1 

This example is interesting because it involves a non-rational dependency on y{x) - the dependent variable 
- thus being out of the scope of most of the symmetry analysis software presently available. It is also 
curious that there are no examples of this type in all of Kamke's set of non-linear second order ODEs. On 
the other hand, using the algorithm here presented, proceeding as in Case A, we determine T, w(x,y), 
and H(y + p(x)) as 

y'(xy' + l)e v 

T = , 2 , / t; w{x,y)=e v : H=l (79) 

y x z + y — 1 

From the last equation we know that we are in Case E. We then determine A and & as in Eqs. (|73) and 
& 

10 We can see this by assuming Ay, = and that Eq. (fn|) does not contain p', and then arriving at a contradiction as 
follows. We first set the coefficients of p' in Eq.(fF)|) to zero, arriving at 



Eliminating ^y'y' gives 



= Bi A + W^y =24y + A x + 9 yly ,y' (A) 



2* y , = B 1 Ay' ~A X 



Differentiating the expression above w.r.t y' and since A y , = we have, 



™ y >y>=B X A 



Finally using Eq.(A), = A, contradicting Fy, ^ 0. 
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A = y'x 2 +y'-l 
y'(xy' + 1) 

* = -2 (80) 

Now, we build Eq.@: 

+ 1 V V ^ + 1 / V V J J B 2 eP 



and, differentiating w.r.t. y', we obtain (Eq.(|77|)), 



2^ + 1- (x 3 +x)y' 2 ( „ l2 \ 2y> -l-2x + xy' r , 

[P +P JH ; — , —g [xp-l) = (82) 



y'^xy' + l) 2 V / (xy' + l) 

Proceeding as explained, dividing by A / and differentiating w.r.t. y' gives 

° L* 2^-1-2* + ^ r ) (xp> - 1) = (83) 

dy' \ {xy' + l)(2xy' + l-(x 3 +x)y' 2 )J 

Solving for p'(x) gives p'(x) — 1/x, from which (Eq.(|46|)): 

V xj y'{xy' + 1) 

Case F 

The final branch occurs when Eq.(fnJ) divided by Kyi does not depend on y' (so that we will not be 
able to differentiate w.r.t y'). In this case we can build a linear algebraic equation for p'(x) as follows. 
Let us introduce the label (3(x,p',p") for Eq.(f77|) divided by Kyi, so that Eq. (f77|) becomes: 

K y ,(x,y') (3(x,p',p") = (85) 

Again, since we obtained Eq.(|77|) by differentiating Eq.(|76|) with respect to y', we see that Eq.(|76|) can 
be written in terms of j3 by integrating Eq.([S5|) with respect to y': 

A(x : y>nx : pi,p") +1 (x,p',p") = B Jl )Bl (86) 

where j(x,p' ,p") is the constant of integration, and can be determined explicitly in terms of x, p' and 
p" by comparing Eq.(86) with Eq.(|7^). Taking into account that (3(x,p',p") — 0, we see that Eq.([S6|) 
reduces to: 

^P'^ = ^§kB7 ( 87 ) 

We can remove the unknowns Bi and B3 after multiplying Eq.(|87|) by e p ^ Bl , differentiating with respect 
to x, and then dividing once again by e p ( x ' Bl . We now have our second equation for p' , which we can 
build explicitly in terms of p', since we know j(x,p',p") and B±: 

^L+B lP ' 1 = 
ax 



Eliminating the derivatives of p 1 between Eq.(|S5|) and Eq.(J88|) leads to a linear algebraic equation in p'. 
Once we have p', the determination of T{x,y') follows directly from Eq.(W6). □ 
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3.2 Integrating factors of the form fj,(y, y') 



Just as in the previous section, from Eq.(16), the ODE family admitting an integrating factor of the form 



y" 



Ky 



For this ODE family, it would be possible to build a pattern matching routine as done in the previous 
section for the case /x(ar, y'). However, it is straightforward to notice that under the transformation 
y(x) — ► x, x — » y(x), Eq.(|89|) transforms into an ODE of the form Eq.([l9]) with integrating factor 
/j,(x, y' ) jy' . This means that the above pattern can be matched by merely changing variables in the 
given ODE and matching Eq.@. It follows that any explicit 2nd order ODE having an integrating 
factor of the form fj,(y, y') can be reduced to a first order ODE by first changing variables, and then using 
the scheme outlined in the previous section (unless the resulting ODE is linear). 



Example: 



/2 

y" - — + sin(x) y' y + cos(x) y 2 = (90) 

y 

Changing variables as in y{x) — > x, x — > y(x) we obtain 

y" H sin(y) y' x — cos(y) x 2 y' = (91) 

x 

Using the algorithm outlined in the previous section, an integrating factor of the form [i(x, y') for Eq.(|9~l|) 
is given by 



y x 



(92) 



from where an integrating factor of the form /i(y, y') for Eq.(|90|) is 1/y, leading to the first integral 

sin(x)y + — + d = 0, (93) 

y 

which is a first order ODE of Bernoulli type. The solution to Eq.(p0|) then follows directly. This example 
is particularly interesting since from |)| we know ODE Eq.(p0|) has no point symmetries. 

3.3 Integrating factors of the form fi(x,y) 

For completeness, we review here the determination of integrating factors of the form fJ,(x, y) for second 
order ODEs, already found in the literature (see for instance Lemma 3.8 in (l(J). Contrary to the cases 
/j,(x, y') or n(y, y'), an integrating factor depending only on x and y can easily be found - when it exists 
- by directly solving the determining equations (|l^) and (pr|). 

^From Eq.(fl6|), the general second order ODE having an integrating factor jj,(x, y) takes the form 

y" = a(x, y) y' 2 + b(x, y) y' + c{x, y), (94) 

where 

a(x,y) = -^, b {x tV ) = -9v±J±, c{x , y)= -^ (95) 
M M M 

for some unknown function G(x,y). As a shortcut to solving Eqs. ( |To| ) and (|Tl]), one can directly tackle 
Eqs.(p5|); the calculations are straightforward. There are two cases to be considered. 

Case A : 2a x — b y ^ 

Defining the two auxiliary quantities 

tp = c y - ac- b x , T = a x . x + a x b + ip y (96) 
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an integrating factor of the form (i(x, y) exists only when 



T 



and is then given by 



fj,(x, y) = exp 



Ox 



bT 



ady I dx — / ady 



(97) 



(98) 



Case B : 2a x — b 
Redefining ip 



y u 

c y — ac, 



an integrating factor of the form (i(x, y) exists only when 
a x ,x - a x b - ip y = 0, 



Then, /x(x, y) is given by 



a dy 



fj,(x,y) = v{x) e 

where v(x) is either one of the independent solutions of the second order linear ODE 



v" = A(x) v' + B{x) v, 



and 



1 = ~dx i ady ' 



A{x) =21-6, 



B(x) = <p+[I- — ) (6-2) 



(99) 



(100) 



(101) 



(102) 



It should be noted that when the attempt to solve the linear ODE Eq.( lOl ) is successful, using each of 
its two independent solutions for integrating factors leads to the general solu tion of Eq.(94), instead of 
just a reduction of order. Also, when the original ODE was linear, Eq.(lOl) is just the corresponding 
adjoint equation, as was to be expected (see for instance Murphy's book). 

3.4 The Connection to PDEs 

Let R{x, y, y') be a first integral of Eq. (|2tj|) . We rewrite Eq.(|l3|) by renaming y' = z 



OR dR _ ,dR „ 

ir + z ^r + V' z ^ir = 

ox ay oz 



(103) 



From Theorem 3.1, if a given PDE of the form Eq.( |103| ) has a particular solution of the form R(x,y, z) = 
F(x, z) + G(x, y), such that R{x, y, z) is nonlinear in y or z; or R(x, y, z) = F(y, z) + G(x, y), such that 
R{x,y, z) is nonlinear in y or z^ 1 , then F and G can be determined in a systematic manner. 

Although this is a natural consequence of the previous sections, it is worth mentioning that the 
dete rmin ation of R using the scheme here presented does not require solving the characteristic strip of 
Eq.( |l03[) , thus being a genuine alternative. 



4 Integrating factors and symmetries 

The main result being presented in this paper is a systematic algorithm for the determination of inte- 
grating factors of the form /j,(x,y') and fi(y,y') without solving any auxiliary differential equations, and 
this last fact is the most relevant point. Nonetheless, it is interesting to briefly review the similarities 
and differences between the standard integrating factor (fj,) and symmetry approaches, so as to have an 
insight of how complementary these methods can be in practice. 

To start with, both methods tackle an n th order ODE by looking for solutions to a linear n th order 
determining PDE in n + 1 variables (see sec. ^|). Any given ODE has infinitely many integrating factors 
and symmetries. When many solutions to these determining PDEs are found, both approaches can, in 
principle, give a multiple reduction of order. 

In the case of integrating factors there is one unknown function, while for symmetries there is a pair 
of infinitesimals to be found. But symmetries are defined up to an arbitrary function, so that we can 
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always take one of these infinitesimals equal to zero; hence we are facing approaches of equivalent levels 
of difficulty and actually of equivalent solving power too. 

Also valid for both approaches is the fact that, unless some restrictions are introduced on the functional 
dependence of fi or the infinitesimals, there is no hope that the corresponding determining PDEs will 
be easier to solve than the original ODE. In the case of symmetries, it is usual to restrict the problem 
to ODEs having point symmetries, that is, to consider infinitesimals depending only on x and y. The 
restriction to the integrating factors here discussed is similar: we considered /i's depending on only two 
variables. 

At this point it can be seen that the two approaches are complementary: the determining PDEs 
for \l and for the symmetries are different^, so that even using identical restrictions on the functional 
dependence of \x and the infinitesimals, problems which may be untractable using one approach may be 
easy or even trivial using the other one. 

As an example of this, consider Kamke's ODE 6.37, appearing in this paper as Eq. (f?8|) : 

y" + 2yy' + f(x) (y' + y 2 ) - g{x) = 

As mentioned in the exposition, for arbitrary f(x) and g(x), this ODE has an integrating factor depending 
only on x, easily determined using the algorithm presented. Now, for non-constant f{x) and g(x), this 
ODE has no point symmetries, that is, no solutions of the form y), rj(x, y)], except for the particular 
case in which g(x) can be expressed in terms of f(x) as inp| 



fii o f fi f3 



C 2 exp ^-3/2 J f{x)dx 
4 ^2C*i + J exp ^-1/2 J f{x)dx^j dx 



Furthermore, this ODE has no non-trivial symmetries of the form [£{x,y'), r](x,y')] either, and for 
symmetries of the form [£(?/, y'), rj(y, y')] the determining PDE does not even split into a system. 
Another ODE example of this type is found in a paper by || (1988): 

y"- — -g(x) Py Py'~g'y p+1 = (105) 
y 

In that work it is shown that for constant p, the ODE above only has point symmetries for very restricted 
forms of g{x). For instance, Eq. (|90|) is a particular case of the ODE above and has no point symmetries. 
Nonetheless, for arbitrary g(x), Eq.( |l05| ) has an obvious integrating factor depending on only one variable: 
1/y, leading to a first integral of Bernoulli type: 

^-g(x) 1 f + C 1 = (106) 

y 



so that the whole family Eq.(105) is intcgrable by quadratures 



having integrating factors of the forms n(x) 



We note that Eq.(|28|) and Eq.(105) are respectively particular cases of the general reducible ODEs 



y" = -^±^y'--^ (107) 



where fi{x) and G(x,y) are arbitrary; and /i(y): 



11 We are considering here ODEs of order greater than one. 

12 To determine g(x) in terms of f(x) we used the standard form Maple package by Reid and Wittkopf complemented 
with some basic calculations. 

13 To obtain the geiieral ODE family reducible by a given integrating factor we used the routine redode also presented 



in this paper in sec. 
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In turn, these are very simple cases if compared with the general ODE families Eq.(|19|) and Eq. (jS9|) , 
respectively having integrating factors of the forms (i(x, y') and n(y, y'), and which can be systematically 
reduced in order using the algorithm here presented. 

It is then natural to conclude that the integrating factor and the symmetry approaches can be useful 
for solving different types of ODEs, and can be viewed as equivalently powerful and general, and in 
practice complementary. Moreover, if for a given ODE, an integrating factor and a symmetry are known, 
in principle one can combine this information to build two first integrals and reduce the order by two at 
once (|| , chap. 3). 



5 Tests 

After plugging the reducible-ODE scheme here presented into ODEtools, we tested the scheme and 
routines using Kamke's non-linear 246 second order ODE example^} The purpose was to confirm the 
correctness of the returned results and to determine which of these ODEs have integrating factors of the 
form fj,(x, y') or fj,(y, y'). The test consisted of determining fx and testing the exactness condition Eq.(^) 
of the product \x times ODE. 

We then ran a comparison of performances in solving a related subset of Kamke's examples using dif- 
ferent computer algebra ODE-solvers (Maple, Mathematica, MuPAD and the Reduce package Convode). 
The idea was to situate the new scheme in the framework of a sample of relevant packages presently 
available. As a secondary goal, we were also interested in comparing the solving performance of the new 
scheme with the one of the symmetry scheme implemented in ODEtools. 

Finally we considered the table of integrating factors for second order non-linear ODEs found in 
Murphy's book and the answers for them returned by all these ODE-solvers. 



5.1 The reducible-ODE solving scheme and Kamke's ODEs 

To run the test with Kamke's ODEs, the first step was to classify these ODEs into: missing x, missing 
y, exact and reducible, where the latter refers to the new scheme. The reason for such a classification is 
that ODEs missing variables are straightforwardly reducible, so they are not the relevant target of the 
new scheme. Also, ODEs already in exact form can be easily reduced after performing a simple check 
for exactness; before running the tests all these ODEs were rewritten in explicit form by isolating y". 
For classifying the ODEs we used the odeadvisor command from ODEtools. All the integrating factors 
found satisfied the exactness condition Eq.(Q). The classification we obtained for these 246 ODEs is as 
follows 



Classification 


ODE numbers as in Kamke's book 


99 ODEs are missing x 
or missing y 


1, 2, 4, 7, 10, 12, 14, 17, 21, 22, 23, 24, 25, 26, 28, 30, 31, 32, 40, 42, 43, 45, 
46, 47, 48, 49, 50, 54, 56, 60, 61, 62, 63, 64, 65, 67, 71, 72, 81, 89, 104, 107, 
109, 110, 111, 113, 117, 118, 119, 120, 124, 125, 126, 127, 128, 130, 132, 
137, 138, 140, 141, 143, 146, 150, 151, 153, 154, 155, 157, 158, 159, 160, 
162, 163, 164, 165, 168, 188, 191, 192, 197, 200, 201, 202, 209, 210, 213, 
214, 218, 220, 222, 223, 224, 232, 234, 236, 237, 243, 246 


13 arc in exact form 


36, 42, 78, 107, 108, 109, 133, 169, 170, 178, 226, 231, 235 


40 ODEs are reducible 
with integrating factor 
fJt(x,y') or n{y,y') and 
missing x or y 


1, 2, 4, 7, 10, 12, 14, 17, 40, 42, 50, 56, 64, 65, 81, 89, 104, 107, 109, 110, 
111, 125, 126, 137, 138, 150, 154, 155, 157, 164, 168, 188, 191, 192, 209, 
210, 214, 218, 220, 222, 236 


28 ODEs are reducible 
and not missing x or y 


36, 37, 51, 66, 78, 97, 108, 123, 133, 134, 135, 136, 166, 169, 173, 174, 175, 
176, 178, 179, 193, 196, 203, 204, 206, 215, 226, 235 



Table 1. Missing variables, exact and reducible Kamke's 246 second order non-linear ODEs. 



^From the table above, ~ 30% of these 246 ODEs from Kamke's book are reducible to first order 
using the scheme here presented. Also, although the symmetry scheme implemented in ODEtools - which 
works with dynamical symmetries and includes heuristic procedures - finds symmetries for 191 of these 
246 ODEs, it is unsuccessful in reducing the order of five ODEs which the new scheme does reduce. 
These are the ODEs numbered 36, 37, 123, 215, and 235. It is interesting to note that ODEs 36, 37 and 
123 have no point symmetries; ODE 215 lead to a third order PDE system whose solution - in terms of 

14 Kamke's ODEs 6.247 to 6.249 cannot be made explicit and are then excluded from the tests. 
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elliptic integrals - can be obtained by computer plus hand if one uses trial and error; and for ODE 235, 
the determining PDE for the symmetries does not split into a system due to the presence of an arbitrary 
function of y'. Also, none of the other computer algebra ODE-solvers tested during this work succeeded 
in solving or reducing the order of any of the se fi ve ODEsp|, even though the corresponding integrating 



factors depend on only one variable (see sec. 5.2) 



For ODE 215, which we write in explicit form as 



ode: = y"= : 3 i \ -JW (109) 



,,, _ (6?/-f)?/ 2 
4 y 3 — a y — b 

the integrating factor found by the computer algebra implementation of the new scheme (see sec. § is0: 

> mu = intf actor(ode,y(x)) ; 

M =i (110) 

y 

This integrating factor leads to a reduced ODE which can be solved as well, resulting in the following 
implicit solution in terms of an elliptic integral^: 

> odsolve (ode) ; 

f(x)dx 

dx-Cx ^dz + C 2 = (111) 

J V-4z 3 + za + b 



The integrating factor found for ODE 37 (see Eqs.(p8|) and (|32|)) leads to a reduction of order resulting 
in the most general first order Riccati type ODE; in this example odsolve just returns the reduction 
of order obtained by the new integrating factor scheme. For ODE 123, the integrating factor found is 
1/y, and the reduced ODE is also a generic Riccati ODE. Finally, ODE 235 appears in Kamke's book 
written in exact form, but the ODE is interesting because it contains three arbitrary functions, of the 
first derivative, the dependent and the independent variables, respectively. Such an arbitrary dependence 
makes this ODE almost intractable for most computer algebra ODE-solvers and related packages. We 
then first isolated the highest derivative as to make the ODE non-exact 



odc:= ^ jw^i (112) 



The integrating factor here found is 
> mu = intf actor (ode ,y(x) ) ; 



fi = H(y') (113) 



Concerning timings, it is worth mentioning that in the specific subset of 28 Kamke's examples which 
are not missing variables, the average time consumed by odsolve in solving each ODE using the new 
scheme was 2.5 sec, while using symmetries this time jumps to 21 sec. These tests were performed using a 
Pentium 200, 64 Mb RAM, running Windows 95. In summary: for these 28 ODEs having an integrating 
factor of the form /j,(x, y') or /j,(y, y'), the new scheme seems to be, on the average, ~ 10 times faster than 
the symmetry scheme. 



5 A table of results obtained using the Reduce package CRACK, kindly sent to us by Dr. T. Wolf, shows that out of 
these five, CRACK is determining symmetries only for ODE 215, and concluding that there exist no point symmetries for 
ODEs 6.36, 6.37, 6.123 and 6.235, but perhaps for some undetermined special values of the function parameters entering 
the ODEs (see details for ODE 6.37 in sec. [§. 

16 In what follows, the input can be recognized by the Maple prompt >. 

17 For ODE 6.215, there is a typographical mistake in Kamke's book concerning the reduced ODE: instead of 
. . . yj 4j/ 3 — #22/ — 33 ••• , one should read ... \/ iy 3 — 32J/ — 33- ■ ■ ■ 
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5.2 Comparison of performances 



With the classification presented in Table 1. in hands, we used different computer algebra systems to 
run a comparison of performances in solving these ODEs having integrating factors of the form /i(x, y') 
or /i(y,y'). For our purposes, the interesting subset is the one comprised of the 28 ODEs not already 
missing variables (see Table 1.). The results we obtained are summarized in the following tableFj: 



Kamkc's ODE numbers 




Convode 


Mathematica 3.0 


MuPAD 1.3 


ODEtools 


Solved: 


51, 166, 173, 174, 
175, 176, 179. 


78, 97, 108, 166, 
169, 173, 174, 175, 
176, 178, 179, 206. 


78, 97, 108, 133, 
166, 169, 173, 174, 
175, 176, 179, 


51, 78, 97, 108, 133, 
134, 135, 136, 166, 169, 
173, 174, 175, 176, 178, 
179, 193, 196, 203, 204, 
206, 215. 


Totals: 


7 


12 


11 


22 


Reduced: 








36, 37, 66, 123, 226, 
235. 


Totals: 
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Table 2. Performances in solving 28 Kamke's ODEs having an integrating factor fi(x,y') or /i(y,y') 

As shown above, while the scheme here presented is finding first integrals in all the 28 ODE examples, 
opening the way to solve 22 of them to the end, the next scores are only 12 and 11 ODEs, respectively 
solved by Mathematica 3.0 and MuPAD 1.3. 

Concerning the six reductions of order returned by odsolve, it must be said that neither MuPAD 
nor Mathematica provides a way to convey them, so that perhaps their ODE-solvers are obtaining first 
integrals for these cases but the routines are giving up when they cannot solve the problem to the end. 

Maple R4 is not present in the table since it is not solving any of these 28 ODEs. This is understandable 
since in R4 the only methods implemented for high order non-linear ODEs are those for ODEs which 
are missing variables. This situation is being resolved in the upcoming Maple R5, where the ODEtools 
routines are included in the Maple library, and the previous ODE-solver has been replaced by odsolvep*[ 

Although the primary goal of this work is just to obtain first integrals for second order ODEs, it is 
also interesting to comment on the six ODEs shown in Table 2. for which the new scheme succeeds in 
determining integrating factors but the reduced ODEs remain unsolved. First of all, for ODEs 36, 37 
and 123 the reduction of order lead to general Riccati type ODEs, so that in these cases no more than a 
reduction of order should be expected. Concerning ODE 235 (Eq.flll2|)), the reduced ODE is: 



V rv r 

H(z)dz+ / G{z)dz+ F{x)dx + Ci=0 (114) 

Methods for solving such a first order ODE are known only for very special explicit combinations of H, 
G and F. Concerning ODEs 66 and 226, the obtained reduced ODEs are the same as those shown in 
Kamke's book, and are out of the scope of odsolve. 



5.3 The reducible- ODE scheme and Murphy's table of integrating factors 

There is an explicit paragraph in Murphy's book concerning integrating factors of the form [i(y ), where 
it is shown a table with four second order non-linear ODE families for which fi(y') is already known. The 
first two families are trivial in the sense that they are already missing variables. The third of these ODE 
families is: 



ode:= y" = P{x)y' + Q{y)y' 2 (115) 

where P and Q are arbitrary functions of its arguments; this is actually Liouville's ODE. The integrating 
factor mentioned in the book is the same found by the scheme here presented: y'\ and the corresponding 
reduced ODE can be solved in implicit form: 

18 When building the statistics for ODEtools, we passed to odsolve the optional argument [reducible] , meaning: try 
the reducible scheme, and if it does not solve the problem just give up. To solve the reduced ODE all of odsolve's methods, 
including symmetries, were used. The inpu t and output in the respective format for all the packages tested arc available in 
littp: / / dft.if. uer j . br / odetools / mu_odcs .zip 

±a However, the scheme here presented was not ready when the development library was closed; the reducible scheme 
implemented in Maple R5 is able to determine, when they exist, integrating factors only of the form n(y'). 



18 



fP(*)dx dx _ / e -fQ(z)dz+C ldz + ( J 2=0 ( U6) 



> odsolve (ode) ; 



The fourth ODE family is the most general second order ODE having 1/y' as integrating factor (see sec. 
I): 

ode:= y" = ^^y' + ^^y' 2 (H7) 

for some function R(x, y). Here the new scheme finds the integrating factor 1/y' and returns the reduced 
ODE 

\n{y')-R{x,y) + C 1 =Q (118) 

actually a generic first order ODE^j. 

6 Computer algebra implementation 

We implemented the scheme for finding integrating factors described in sec. ^ in the framework of the 
ODEtools package j|, taking advantage of its set of programming tool routines specifically designed to 
work with ODEs. The implementation consists of: 

• The plugging of the reducible- ODE solving scheme here presented in the block of methods for nonlinear 

second order ODEs of the ODEtools command odsolve; 

• The extension of the capabilities of the ODEtools intfactor command to determine integrating factors 

for non-linear second order ODEs using the scheme here presented; 

• A new user-level routine, redode, returning the most general explicit ODE having a given integrating 

factor (Eq.(|l|)); 

The computational implementation follows straightforwardly the explanations of sec. ^| and includes 

three main routines, for determining T(x,y'), ^(x) and the reduced ODE R(x, y, y') 7 respectively. Call- 
ings to these routines were in turn added to both the intfactor and odsolve commands, so that the 
scheme becomes available at user-level. 

A test of this implementation in odsolve and some related examples are found in sec. |5|. Since 
detailed descriptions of the ODEtools commands are found in the On-Line help, we have restricted this 
section to a description of the new command redode followed by two examples. 

Description of redode 
Command name: redode 

Feature: returns the reorder ODE having a given integrating factor 
Calling sequence: 

> redode (mu, n, y(x)); 

> redode(mu, n, y(x), R) ; 

Parameters: 

n - indicates the order of the requested ODE. 

mu - an integrating factor depending on x, y, j/^™ -1 ). 

y(x) - the dependent variable. 

R - optional, the expected reduced ODE depending on x, y, j/™ -1 ). 

Synopsis: 

• Given an integrating factor fj,(x,y, ...y( n '), redode's main goal is to return the ODE of order n 
having fi as integrating factor 



20 For the third ODE family, Mathematica 3.0 returns a wrong answer and MuPAD 1.3 gives up, while for the fourth 
family, Mathematica gives up and MuPAD returns an ERROR message. 
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A 4 



d_ 

dx 



where G = G(x, y, ...?/"~ 2 )) is an arbitrary function of its arguments (see sec. ||). This command is 
useful to identify the general ODE problem related to a given fi, as well as to understand the possible 
links between the integrating factor scheme for reducing the order and other reduction schemes (e.g., 
symmetries). 

• When the expected reduced ODE (differential order n— 1), here called R, is also given as argument, 
the routine proceeds as follows. First, a test to see if the requested ODE exists is performed: 

n{x,y,...yt n - r >)=u(x,v,...v<- n -Q) ^-^R(x,y, (119) 

for some function v(x. y, ...j/™ -1 )). If the problem is solvable, redode will then return an reorder 
ODE^ = - $(ar,y, ...,y (n " 1) ) satisfying 

(x(x,y,...,y( n -^) ODE™ = ± (u(x,y, ...y^) R(x, y, y^ 1 ^ (120) 

that is, an ODE having as first integral v R + constant. 

• When the given \x does not depend on y(™ _1 ) and R is non-linear in y( n ~ x \ the requested reorder 
ODE nevertheless exists if R can be solved for z/" -1 ). 

Examples: 

The redode command is interesting mainly as a tool for generating solving schemes for given ODE 
families; we illustrate with two examples. 

1. Consider the family of second order ODEs having as integrating factor fj, = F(x) - an arbitrary 
function - such that the reduced ODE has the same integrating factor. We want to set up an algorithm 
such that, given a second order linear ODE, 

y " = ipi( x )y' + M*)y + Mx) ( 121 ) 

where there are no restrictions on ipi(x), ip2(x) or ip3(x), the scheme determines if the ODE belongs to 
the family just described, and if so it also determines F(x). The knowledge of F(x) will be enough to 
build a closed form solution for the ODE. 

To start with we obtain the first order ODE having F(x) as integrating factor via 

> ode_l := redode(F(x), y(x)); 

o^:=y' = -^{y d -^ + ^l(x)) (122) 

where -Fl(x) is an arbitrary function. To obtain the second order ODE aforementioned we pass ode_l 
as argument (playing the role of the reduced ODE) together with the integrating factor F(x) to obtain 

> ode_2 := redode(F(x), y(x), ode_l) ; 

If, dF(x) d 2 F(x) d_Fl(x)\ , . 

° d62 : = = ~~F(x) ( dx + ^ + ^ ) ^ 

Taking this general ODE pattern as departure point, we setup the required solving scheme by com- 
paring coefficients in Eq.(121) and Eq.( |l23| ), obtaining 

F(x)^r = ^ Fjx)^lx^= Mx) (124) 
By solving the first equation, we get F(x) as 
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^Idx 

F{x) =de J (125) 
and by substituting this result into the second one we get the pattern identifying the ODE family 

5^-3(^)) 2 -^) = o (126) 

Among the ODE-solvers of Maple R4, Mathematica 3.0, MuPAD 1.3 or Convode (Reduce), only those 
of MuPAD and Maple succeed in solving this ODE family. 

2. Consider the second order ODE family having as integrating factor /i = F{x) - an arbitrary function 
- also having the symmetry^] [£ = 0, r\ = F(x)], and such that the reduced ODE is the most general first 
order linear ODE 

odei := y' =A(x)y + B(x) (127) 

where A{x) and B{x) are arbitrary functions. To start with, we obtain the aforementioned second order 
ODE having the integrating factor F(x) as in Example 1. 

> ode_2 := redode (F(x) , y(x), ode_l) ; 



dA{x) 




''" = ' ~dx + ' ''' + ' l( ' ' ^1*' ~^ + ^T <^> 

In this step, ode_2 is in fact the most general second order linear ODE. If we now impose the symmetry 
condition^ A(ode_2)=0, where X = [0, F(x)] we arrive at the following restriction on A(x) 

- F (x) d -f± - 2 ( d Jm A{x) + I^d_ + *m = o (129) 

ax \ ax J b (x) dx z 



Solving this ODE for A(x), introducing the result into Eq.(128) and disregarding the non-homogeneous 



term (irrelevant in the solving scheme) we obtain the homogeneous ODE family pattern: 

y + H(x)y' (130) 



1 dH{x) 3 V dx J 1 —j^- 



2 dx 4 (H(x)f 2 H{x) 

where we introduced H(x) = (F(x))~ 2 . Although this ODE family appears more general than the one 
treated in Example 1., the setting up of a solving scheme here is easier: one just needs to check if the 



coefficient of y in a given ODE is related to the coefficient of y' as in equation Eq.(130), in which case 
the integrating factor is just , 1 

y/H(x) 

7 Conclusions 

This paper presented a systematic method for obtaining integrating factors of the form /i(x, y') and fx(y, y') 
- when they exist - for second order non- linear ODEs, as well as its computer algebra implementation in 
the framework of the ODEtools package. The scheme is new, as far as we know, and the implementation 
has proven to be a valuable tool since it leads to reductions of order for varied ODE examples, as shown 
in sec. |[ Actually, the implementation of the scheme solves ODEs not solved by u sing standard or 
symmetry methods (see sec. |]) or other computer algebra ODE-solvers (see sec. 5.2); furthermore, it 



involves only algebraic operations, so that - in principle - it gives answers remarkably faster than the 
symmetry scheme. 

21 Here we denote the infinitesimal symmetry generator by [£,77] = + 

22 For linear ODEs, symmetries of the form [0,F(x)] are also symmetries of the homogeneous part. 
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It is also worth mentioning that restricting the dependence of /i in Eq.(Q) to fj,(x,y'), does not lead to 
a straightforwardly solvable problem except for few or simple cases. Moreover, when the determination 
of /j, from Eq.(Q) is frustrated, there is no way to determine whether such a solution /i(x,y') exists. It 
is then a pleasant surprise to see such integrating factors - provided they exist - being systematically 
determined in all cases and without solving any differential equations, convincing us of the value of the 
new scheme. 

On the other hand, we are restricting the problem to the universe of second order ODEs having 
integrating factors depending only on two variables - the general case is fj,(x,y,y') - and even so, for 



integrating factors of the form n(x,y) the method may fail in solving the auxiliary linear ODE Eq.(LOl) 
which appears in one of the subcases. 

Some natural extensions of this work then would be to develop a scheme for building integrating 
factors of the forms considered in this work, now for higher order ODEs, at least for restricted ODE 
families yet to be determined. Concerning these extensions, the redode routine presented, designed to 
find the most general n'^order ODE having a given integrating factor, optionally reducing to a given k th 
order ODE [k < nY, can be of use in investigating further problems. We are presently working on these 
possible extensions^ and expect to succeed in obtaining reportable results in the near future. 
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Appendix A 



This appendix contains some additional information which may be useful as a reference for developing 
computer algebra implementations of this work, or for improving the one here presented. 

As explained in sec. 3.1.2 , the scheme presented can be subdivided into six different cases: A, B, C, 
D, E and F. Actually, there are just five cases since case B is always either A or C. From the point of 
view of a computer implementation of the scheme it is interesting to know what one would expect from 
such an implementation concerning Kamke's ODEs and the cases aforementioned. We then display here 
both the integrating factors obtained for the 28 Kamke's ODEs used in the tests (see sec. |J) and the 
"case" corresponding to each ODE. 



Integrating factor 



Kamke's book ODE-number 



Case 



1 



/'/(*)<* 



36 



37 



51, 166, 169, 173, 175, 176, 179, 196, 203, 204, 206, 215 



b+y 



/2\3/2 



66 



78 



97 



108 



123 



133 



y -i 
(i+y')y' 



134 



V-l 



135 



y -i 
h(y') 



136 



2 xy —1 

(i + yT 1 



174 



178 



y' (i+2j/y') 



193 



226 



h(y') 



235 



Table A.l Integrating factors for Kamke's reducible and ODEs not missing variables. 
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